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We analyze the neutrino mass matrix entries and their correlations in a probabilistic fashion, 
constructing probability distribution functions using the latest results from neutrino oscillation fits. 
Two cases are considered: the standard three neutrino scenario as well as the inclusion of a new 
sterile neutrino that potentially explains the reactor and gallium anomalies. We discuss the current 
limits and future perspectives on the mass matrix elements that can be useful for model building. 

PACS numbers: 



The year 2012 represents a milestone in neutrino physics. Thanks to the measurement of the last mixing angle 
of the standard neutrino oscillation scenario, 613, by the reactor experiments Double-CHOOZ [1 j, Daya-Bay[2] and 
RENO [3] , after the first positive evidence from accelerators [4] [5] , the mixing in the leptonic sector is starting to shape 
up. The impact of a rather unexpectedly large mixing angle #13 is two fold: it promotes the discovery of CP violation 
in the neutrino sector to a yet daunting but conceivable task, and at the same time it proves that the description of 
neutrino oscillation data must involve all three Standard Model neutrino flavors. 

Hints of the sensitivity to CP phases are already showing their first signs when we combine accelerator — >• v e with 
reactor v e — » v e data [6] or when we perform global fits [3 [8]. Furthermore, the combined fit of neutrino oscillation 
data shows for the first time a very precise and almost complete determination of the parameters that enter the 
standard neutrino oscillation scheme. In fact, in spite of the unknowns (neutrino mass hierarchy, absolute neutrino 
mass scale, CP phases and the correct octant for #23) all measured parameters are now so well determined that it is 
enough to quote them by giving the best fit value with the 1 a uncertainty. 

However, not all neutrino data can be explained by this standard scenario of three flavor neutrinos. In fact, along the 
years a number of so-called anomalies have crept into the picture. First, the excess of v e events in the — ^ v e mode 
observed by the short baseline LSND [9] experiment, now also supported by MiniBOONE data [TP] , gave rise to the 
long-standing LSND anomaly. Second, the deficit of v e compared to expectations observed by the source calibration 
experiments performed in the gallium radiochemical solar neutrino detectors GALLEX [11] and SAGE [12]. This is 
the so-called gallium anomaly. Third, and more recently, a re-evaluation of the reactor v e flux [13] [14] resulted in 
an increase of the total flux by 3.5%. While this increase has essentially no impact on the results of long baseline 
experiments, it induces a deficit of about 5.7% in the observed event rates for short baseline (< 100 m) reactor 
neutrino experiments. This problem has been referred to as the reactor antineutrino anomaly [T5] , 

There are attempts in the literature that try to explain some or all of these anomalies by extending the standard 
picture to include one or more sterile neutrinos [T6HT8] . These extensions, as a rule, cannot make appearance and 
disappearance experiments compatible. However, if one disregards the anomaly connected with the appearance 
experiments LSND/MiniBOONE (for instance, assuming it is not due to oscillations), it is possible to construct a 
coherent picture of all solar, atmospheric, reactor and accelerator neutrino oscillation data adding one extra sterile 
neutrino to the standard framework. This constitutes what has been known as the 3+1 scenario. 

Given the current status of the mixing parameters measurements, and in view of the progress expected in the near 
future, we think it is timely to analyze the possible structures and correlations among the neutrino mass matrix 
elements that are compatible with data. In this sense, we update Refs. [19] [20] using the most recent available data. 
However, our analysis will be probabilistic, since we will construct probability distribution functions for each element 
of the neutrino mass matrix. We also discuss how this might change with better determination of the presently 
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known oscillation parameters, as well as, the Dirac CP-phase 5. We hope this can be helpful to understand better 
the patterns statistically preferred by data, and serve as a guide for model builders. 

We organize our paper as follows. In Sec. [IT] we describe how we will proceed for the construction of probability 
density functions (PDF) for each element of the neutrino mass matrix and what are the assumptions in each case. 
In Sec. Ill we analyze the possible textures of the mass matrix in the standard scenario, discussing the correlations 
among matrix elements in the hierarchical and almost degenerate cases. We also discuss the future prospects for 
better determining these matrix elements with neutrino oscillation and non-oscillation data, and the possible impact 
on the theory. In Sec. IV we extend our analysis to include the possibility of a sterile neutrino with mass and mixings 
allowed by the reactor and gallium anomalies. We discuss what is the mass matrix pattern in this scenario and how 
different it will be from the standard case. Finally, in SecjVl we make our last comments and draw our conclusions. 



II. RECONSTRUCTING THE NEUTRINO MASS MATRIX FROM EXPERIMENTAL DATA 

To access the impact of the progress on the determination of the neutrino oscillation parameters in the last year on 
the knowledge of the low energy effective neutrino mass matrix, in a probabilistic way, we will construct a PDF for 
each element of the mass matrix in the gauge basis, 



m aP = Y j m i e-^U* ai Ul i (1) 

i 

with a,/3 = e, /i, r, U a i being the elements of the mixing matrix, mi the neutrino masses and A^ the Major ana- type 
CP phases, using the most recent available information from the combination of neutrino oscillation data. Without 
loss of generality we will take A2 = in our parametrization. 

We will use the following best fit points for the standard mixing parameters [8J 

(7.50 ±0.185) x HT 5 eV 2 
0.30 ±0.013 

(±2.47 ± 0.07) x 10" 3 eV 2 (normal ordering) 
(-2.43 ± 0.06) x 10" 3 eV 2 (inverted ordering) 
0.023 ± 0.0023 
f 0.41 ±0.037 (1 st octant) 
\ 0.59 ± 0.022 (2 nd octant) 

and assume these parameters to follow a normal distribution with mean at the best fit point and standard deviation 
equal to the 1 a uncertainty. We will take the unknown Dirac-type and Majorana-type CP phases to be flat distributed 
between and 2tt. 

For the 3+1 scenario we will fix the squared mass difference between the sterile and the lightest state to two different 
experimentally allowed values, Am^ = 1.71 eV 2 and Am 2 2 = 0.95 eV 2 , allowing the corresponding mixing to vary 
with a flat distribution inside the ranges \U e 4\ 2 = 8 x 10 -3 — 4 x 10 -2 and |Z7 e 4| 2 = 8 x 10 -3 - 2.5 x 10 -2 , respectively. 

To construct the PDF of each m a p we use a Monte Carlo method. For all the mixing parameters, we generate 
random numbers according to their assumed distribution, and compute all the elements of m a p in each case. In this 
manner their distribution will be naturally correlated. 

Since we do not know the neutrino mass hierarchy, the correct octant for #23 and the absolute neutrino mass scale, 
we will have to analyze each case separately. 



sin 2 #12 = 

Am 2 31 = 
Am 2 2 

sin 2 #13 = 

sin 2 #23 



III. THE STANDARD SCENARIO 

In the standard scenario we use the standard parametrization for the mixing matrix, 

(C12C13 S12C13 s 13 e~ lS \ 

-S12C23 - c 12 s 13 S23e l6 C12 c 23 - S12S13S23 e l6 c 13 s 23 , (2) 
Sl2 S23 ~ C12S13C23 e%b ~ C 12 S 23 - S12S13C23 e lS C13C23 / 

with Sij = sin 6^-, Cij = cos 6ij and S the Dirac-type CP phase. The two additional Majorana-type CP phases are 
denoted by Ai and A3, and the three neutrino mass eigenstates are ordered either as mi < < m 3 (normal ordering) 
or as 777,3 < mi < ^2 (inverted ordering). 
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We will study six different cases: hierarchical with mi — >• and 623 both in the first and second octant; hierarchical 
with 777,3 —> and 023 in the first and second octant; quasi-degenerate with mi ~ 7772 ~ ms ~ 0.1 eV and 023 in the 
first and second octant. There are two possible ordering also in the quasi-degenerate case; however, we have checked 
that the results are very similar. 

In Fig. [l]we illustrate, for the case m\ — > 0, the impact of the determination of sin 2 #13 on the PDF of \m ee \. 
The distribution labeled "before" (magenta) is obtained assuming sin 2 #13 to be flat distributed between and 0.04 
(CHOOZ limit [21 j), while the one labeled "after" (blue) shows the current situation. The two peaks in the "after" 
distribution are due to the interference between the real U^ 2 m 2 term and the complex U^ 3 ms term (see Appendix |a| 
for detailed expressions), which depends on the cosine of the randomly distributed CP-phases. This term depends on 
013 which is not compatible with zero anymore and thus gives a sizable contribution. The distance between the peaks 
depends on 7773: a larger 7773 would place the peaks further apart. 

In Fig. [2] we illustrate, again for the case mi — » 0, the effect of the determination of the 023 octant on the PDF 
of I^VaJ- O n the left (right) panel we show the current distribution of |ra M/x | for the best fit value of sin 2 2 3 in the 
first (second) octant. The asymmetric two peaks structure is due to the fact that #12 is not maximal. A larger 7773 
would shift the right endpoint of the distribution to higher values of |m M/i |, while a larger 7772 would separate the two 
peaks. We observe that the distribution on the right panel is thinner than the one on the left panel, but this is simply 
because the uncertainty on sin 2 023 is smaller in the second octant. For comparison, in both panels we also show the 
PDF taking the MINOS result sin 2 23 = 0.5 ± 0.1 [221 ■ 
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FIG. 1: PDF for \m ee \ when mi —> 0. The "before" ("after") distribution corresponds to the situation before (after) the 
determination of sin 2 613 by the reactor experiments. 




FIG. 2: PDF for \m^\ when mi — )> 0. The left (right) panel corresponds to the distribution for 2 3 in the first (second) octant. 
We also show the distribution for the case of maximal mixing with MINOS uncertainty [22] . 
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in meV 


Element 


mi -+ (1 st oct.) 


mi ^ (2 nd oct.) 


m 3 (1 st oct.) 


m 3 (2 nd oct.) 


\m ee \ 


1.3 - 4.2 


1.4 - 4.2 


18 - 53 


18 - 54 


\m e ^\ 


1.5 - 8.8 


2.5 - 9.4 


2-40 


2-35 


\m eT \ 


2.7-9.3 


1.5 - 9 


2-35 


2-40 


\m^\ 


14 - 27 


24 - 34 


6-35 


2-25 


% T 


21 - 28 


20.5 - 28 


9- 27 


9-27 


\m TT \ 


23 - 35 


14 - 26 


2-26 


6-34 



TABLE I: Range of allowed values of |ra a/ g| at 95.45 % CL for the very hierarchical cases. 

A. Correlations Among Matrix Elements 

1. Hierarchical Case 



In the case mi — > 0, 777,3 ~ 0.05 eV ^ 777,2 ~ 0.009 eV and only two CP-phases, S and A3, are relevant. Due to 
the [I — y t symmetry, accomplished by S23 — » C23 and C23 — » — S23, the PDFs for the solution in the first #23 octant 
are basically the same as for the solution in the second octant (apart from the uncertainty in the determination of 
sin 2 #23, which is smaller for the second octant) as long as we replace: \m er \ |ra 



\m u 



In Fig. m\ we show the correlations among the absolute values of some of the matrix elements m a p for mi — » and 
#23 in the first octant. In appendix [B] we present a complete set of those plots, as well as the corresponding plots for 
#23 taken to be in the second octant. This was done by constructing a two-dimensional PDF for each pair of elements. 
We also present in each case on the top and to the right of each two-dimensional distribution the projected PDFs. 
In these plots we use blue, green and red for the allowed region at 68.27%, 95.45% and 99.73% CL, respectively. The 
range of the values allowed at 95.45% CL are given in Table [l| 



From Figs 
pairs |ra* J x 



1 

: \rr 



is not very correlated to any other element. However, due to #23 the 



x |ra M/i |, \m^\ x |ra Mr |, |m MT | x |ra rr |, \m^\ x |ra rr | and |ra e/i | x \m TT \ are very correlated. 



13 and 14 we observe that \m ee 
W I? |ra er | x Ira^^l, \rn j j LjJL \ x |ra MT |, \m, jLT ( 
So, for instance, if a model predicts |ra e/i | ~ 6 meV, \m eT \ must be in the range 5-8 meV, while without taking into 
account this correlation the allowed range would be 2.7-9.3 meV at 95.45% CL. 

From the expressions for the matrix elements given in appendix [A] it is easy to show that in this case: 



ra, 



12 



e/x| 



ra fi 



2 2 
m 3 y 



2^2 



m 9 x 



|ra rr | 2 ~ra 2 (l-2z 2 ) + | 



and 



ra 



•/z/x| 



\m^ r \ 2 + |ra rr | 2 - ra 2 (l - z 2 ), 

where x = sin#i 2 , y = sin #13 and z = sin #23 • 

Due to the prevalence of the ra3 mass, the determination of sin 2 #13 with the present uncertainty of 10% by the 
reactor experiments not only affected the range of |ra ea |, a = e,/i,r, but also their PDFs. Also, the determination of 
sin 2 #23 with an uncertainty of 9% (4%) in the first (second) octant changes the range and the form of the PDFs of 
all matrix elements except for |ra ee |. 

In the case ra3 — >• 0, mi « ra2 ~ 0.05 eV and only two CP-phases, S and Ai, are important. Again the PDFs 
obtained for the absolute value of the matrix elements in the first and second octant of #23 are basically the same as 
long as we replace fi <H> r. In Fig. [4] we show two-dimensional PDFs for some pairs of elements of the matrix m a p 
in the case ra3 — >• and #23 in the first octant. A more complete set of plots can be found in appendix |Bj while in 
Table [l| we present the 95.45% CL allowed ranges. 

Generally, the predominant terms comprise rai or 777,2, which have similar sizes, and their contributions involve #12 
and #23, which are not maximal, without being suppressed by #13. There are at least three consequences of such a 
fact. First, the determination of sin 2 #13 by the reactor experiments basically did not affect the range of the matrix 
elements but changed the shape of some of their PDFs. Second, the determination of sin 2 #23 with an uncertainty of 



5 




1 2 3 4 5 10 20 30 20 25 30 

|m ee | (meV) \m m \ (meV) |m^-| (meV) 



FIG. 3: PDFs for the distribution of the absolute value of several pairs of matrix elements. Top panels: \m ee \ x \m^ T \ (left), 
\fTien\ x \m eT \ (center) and Ira^l x \m eT \ (right). Bottom panels: \m ee \ x \m e ^\ (left), |m M/x | x \m TT \ (center) and \m^ T \ x \m TT \ 
(right). At the top and right of each two dimensional PDF we show the PDF of the absolute value of the corresponding matrix 
element. In these plots we use blue, green and red for the allowed region at 68.27%, 95.45% and 99.73% CL, respectively. Here 
mi —> and 623 is assumed to be in the first octant. 



9% (4%) in the first (second) octant changes the range |m e/x |, |m er |, \m^\ and |m rr |, while the shape of the PDFs 
remain basically the same except in the case |m M/x | and \m TT \. Last, the mass matrix entries are very correlated. 
We observe the strong correlations between all pairs of elements in Figs. [4j [15] and 16 So in this case it is even 



more important to take into account these correlations in model building. For instance, we can easily see, from the 
expressions in appendix [Aj that in this case 

\m^ r \ ~ y/z 2 (l- z 2 ) |m ee |, 
\m ee \ 2 ~ ml - (1 - z 2 ) _1 |m e/i | 2 , 

and 



i 1 A 1 -* 2 ) . 1 

\m eil \ ~ \ - 2 — \m er \. 



This behavior is confirmed by Fig. [4j 



2. Quasi- Degenerate Case 

In the quasi-degenerate case, mi ~ rri2 ~ and the effect of the ordering is very small. In this case, all masses 
and CP phases play a role. 

As an example, we take mi = 0.1 eV. In Fig. [5] we show the correlations among the absolute values of some of 
the matrix elements for the normal mass ordering and #23 in the first octant. In appendix [B] one can find a 
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FIG. 4: Same as Fig. [3] but for rriz — >> and #23 in the first octant. 





in 


meV 






Element 


-j^st 


oct. 


2 nd 


oct. 


\m ee \ 


35 


- 108 


36 - 


107 




5 


- 90 


5 - 


80 


m eT 


5 


- 80 


5 - 


90 


\™>w\ 


10 


- 115 


22 - 


-110 


hp 


10 


- 115 


9 - 


110 


\m TT \ 


20 


- 115 


12 - 


115 



TABLE II: Range of allowed values of \m a /3\ at 95.45 % CL for the quasi-degenerate case with mi = 0.1 eV. 



complete set of plots for this case, Fig. 17 , and for #23 in the second octant, Fig. 18 We use the same color coding 
as in previous figures. The range of the values allowed at 95.45% CL are given in Table |TT} 

The correlations here are either similar to the very hierarchical case in normal ordering or in the inverted ordering. 



For example, the PDFs for \m ee \ x |ra e/i |, \m^\ x |m rr |, |m M/x | x \m er \ and |ra e/i | x \m e 
inverted ordering, while \m^ T \ x \m TT \ and \m ee \ x |m Mr | are more like the normal ordering. 



are correlated like in the 



3. Future Perspectives 



To evaluate the effect of a better determination of the mixing parameters on our knowledge of the mass matrix we 
have studied the effect of reducing by half the uncertainty of a parameter at a time, keeping the other parameters at 
their current uncertainties. 

In Fig. [S] we illustrate the effect of a better determination of sin 2 #13 (top left panel), Am^ (top right panel), 
sin 2 621 (bottom left panel) and sin 2 #23 (bottom right panel). We observe that the effect of a better determination 
of any of these parameters is very small. We have verified that this is true for both mass hierarchies and #23 octants. 
The biggest effect comes from a better determination of sin 2 #23, as one could have guessed, but still this only reduces 



7 




|m ee | (meV) \m m \ (meV) \m^\ (meV) 

FIG. 5: Same as Fig.[3]but for mi = 0.1 eV and #23 in the first octant. 

significantly the 3 a region. We do not show the effect of a better determination of Am^ because it is even smaller 
than for the other parameters. 

On the other hand, a measurement of S with an uncertainty of 20°, that can be envisaged according to Ref. [28] for 
long baseline neutrino oscillation experiments, could be significant. To illustrate this we show in Figs. [7][8] the effect 
of a determination of 5 with an uncertainty of 20° for some of the correlations between pairs of mass matrix elements. 
On the top (bottom) panels of Fig. [7] we can see this for \m ee \ x \m efJ ,\ (\m^^\ x|m er |) in the normal ordering for 5 = 0° 
(left), 180° (center) and 270° (rightj^ while on the top (bottom) panels of Fig.[8]we show \m^ T \ x \m TT \ (\m^^\ x |m er |) 
in the inverted ordering for 5 = 0° (left), 90° (center) and 180° (right), These cases for normal (inverted) ordering 
and S = 90° (S = 270°) are not shown because they are very similar to S = 270° (S = 90°). 

For the normal ordering, with mi — »• 0, the determination of 5 will play an important role in the correlation of 
\m ee \ with all other mass matrix elements, but will be more significant for \m efl \ or |m er |. This is due to the fact that 
for these mass matrix elements the leading phase term is the one that accompanies cos [2(5 + A3)], just as for \m ee \. 
This will also affect the correlations involving |m^| and |m rr |, since for them the leading phase terms are, in order of 
importance, the ones that go with cos(2As) and cos(5 + 2A3). However, the correlations with \m^ T \ will only slightly 
change because the leading phase term for this element does not depend on 5. 

For the inverted ordering, with 777,3 —> 0? the determination of 5 will play a bigger role in the PDFs of \m TT \ and 
|m M/x |. This is because, as we can see from their expressions in appendix [A} the leading coefficients of cos(2Ai), 
cos(£± 2Ai) and cos 5 are all of the same order. The PDFs of \m efJL \ and |m er Jwill also be affected because the terms 
that depend on 5 are not negligible in comparison to the leading term that depends on cos(2Ai), however their relative 
importance will depend on the #23 octant. The PDFs for \m ee \ and |m MT | are basically independent of 5, the first 
because ms — » 0, the second because these terms are suppressed by sin 2 #13 or factors of this order. 

We also have checked that the effect of the determination of 5 for the quasi-degenerate case with mi = 0.1 eV is 
smaller than for the hierarchical cases because there are more phases involved. 

In the future we also expect to have information from three different sources: neutrinoless double beta decay 
experiments, beta decay experiments and cosmology. In Fig. [9] we show the current allowed region for \m ee \ as a 
function of the effective electron neutrino mass, mp and of the sum of the neutrino masses [38], J2 m i a ^ 99% CL. 
The region allowed by the normal (inverted) mass ordering is in blue (red) and the recent limit on \m ee \ given by 
KamLAND-Zen [29], \m ee \ < (120 — 250) meV, is shown in gray. Cosmology excludes the magenta region ^rrii > 
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(0.2 — 0.6) eV [31 . Notice that the allowed regions were built from the pdfs constructed from data (except for the 
CP phases, which we assumed to be flat distributed). This is why the m ee —> region is absent, as it is very unlikely 
to have the necessary degree of cancellations. 

The forecast sensitivity on \m ee \ of the most ambitious neutrinoless double beta decay experiments, after 5 years 
of exposure, is 29-73 meV (GERDA phase-3) and 18-39 meV (CUORE) [30]. 

The KArlsruhe TRitium Neutrino mass experiment (KATRIN) will have a sensitivity on the electron neutrino 
effective mass, mp = VEi \Uei\ 2 m\, of 0.2 eV at 95% CL [32]. 

Cosmological limits today still allow for quasi-degenerate neutrino masses; however, this possibility will soon be 
confirmed or ruled out. According to Ref. [31 , many cosmological probes, with different systematics, will reach a 
sensitivity on ra$ of 0.1 eV or lower. Lyman a forest can reach 0.1 eV, lensing of Cosmic Microwave Background 
0.2-0.05 eV, lensing of galaxies 0.07 eV, observations of the redshifted 21 cm neutral hydrogen line 0.1-0.006 eV, 
galaxy clusters and galaxy distribution surveys 0.1 eV. 

If cosmology will point to a quasi-degenerate neutrino spectrum: 

1. the determination of S by future experiments will not modify much the current correlation among mass matrix 
elements; 

2. \m ee \ should be in the reach of most proposed neutrinoless double beta decay experiments, if neutrinos are of 
Majorana nature, ergo the non-observability of the OvPfi would point to Dirac neutrinos; 

3. mp may be in the reach of KATRIN; 

4. the ordering of neutrino masses can be settled by future neutrino oscillation experiments but will not have great 
impact on the neutrino mass matrix; 

5. if \m ee \ is measured we will also be able to constrain much more \m efl \ and |m er |. However, it will be very 
difficult to say something about the Majorana-CP phases. 

On the other hand, cosmology can also place a limit on such that we will know if we have normal ordering 

with mi — » 0. In this case: 

1. \m ee \ will be out of the reach of the proposed 0^/3/3 experiments; 

2. mp will be out of the reach of KATRIN; 

3. the experimental determination of 5 will increase the correlation among mass matrix elements and help to 
determine the structure of the neutrino mass matrix; 

4. the ordering of neutrino masses should be confirmed by future neutrino oscillation experiments. 

It may be the case that cosmology will rule out a quasi-degenerate spectrum, but not the inverted ordering. If this 
happens: 

1. \m ee \ may be in the reach of the proposed Ovfifi experiments; 

2. mp will be out of the reach of KATRIN; 

3. the ordering of neutrino masses should be determined by future neutrino oscillation experiments; 

4. the experimental determination of 5 will increase the correlation among mass matrix elements and help to 
determine the structure of the neutrino mass matrix, specially if the mass ordering is known; 

5. we may be able to say something about one of the Majorana CP phases, if \m ee \ is measured. 

These future advances may thus provide new clues for the understanding of the flavor problem in the lepton sector. 
Some models of neutrino mixing based on discrete flavor groups have predictions that can be tested in the future. For 
instance, the Lin model [23], where the A4 symmetry is broken by additional Z n parities, predicts 

sin 2 #23 = — I — \= sin #13 cos S. 
2 y2 

Another example is the SUSY model based on the flavor symmetry group S4 x Z 4 x U(l) discussed in [24] where 
the relation 

sin 2 O12 = — h sin #13 cos S + (9(sin 2 #13) 
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appears. Both relations can, in principle, be experimentally tested and, if true, they will impose new correlations 
among the neutrino mass matrix elements. 

Other models, such as the one presented in Ref. [25 , are even more predictive. This model, which is based on 
a type-I seesaw framework with an underlying A4 flavor symmetry, can, given a set of vacuum expectation value 
alignments for the flavon fields which break the symmetry, predict neutrino masses, the mass hierarchy, #23 and 5. 

There are also model-independent approaches to the flavor problem in the neutrino sector. In Ref. [26j [27] , relations 
among the mixing parameters were obtained, in the context of discrete flavor symmetries, under general assumptions 
that the flavor symmetry group is of the von Dyck type. Again these relations can, in principle, be experimentally 
tested, and if ratified by experiment induce more correlations among the mixing matrix entries. 



IV. THE 3+1 SCENARIO 



Whether or not one deems this to be a plausible scenario, we still believe it is important to examine what are its 
consequences to the possible textures of the neutrino mass matrix. 
In this case the mixing matrix can be parametrized as 
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Cl2 ^13 C14 


Sl2 Cis C14 


sis c 1A e~ %5 


s 14 


\ 




~S\2 C23 - C12S13S23 e%b 


Cl2 C 2 3 - 5i25i 3 523 e z6 
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S\2 S23 ~ C12S13C23 e%5 


-C12S23 ~ 5 12 5 13C23 e tS 
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-Cl2 C13 Sl4 


-Sl2 Cis «14 


-sissi4e~ %s 


C14 


/ 



(3) 



where we use the same notation as in Eq. (2|. This expression can be readily derived from [34] once we identify S = #14 
and we assume n = (1, 0, 0), i.e. the sterile state mixes only with the electron neutrino. With this assumptions, there 
are no extra Dirac CP phases. Nevertheless, there is one extra Majorana phase, A4. The expression of all the squared 
matrix elements, in the limit c\s ~ C14 — >• 1, can be found in Appendix [A] 

For simplicity, we examine here two cases: (a) Am^ = 1.71 eV 2 and sin 2 #14 = (0.8 — 4.2) x 10 -2 , (b) Am^ = 0.95 
eV 2 and sin 2 #14 = (0.8 — 2.5) x 10 -2 , which are two possible solutions to the reactor and gallium anomalies [T5] . 
These solutions seem at first glance to be at odds with cosmology, but we will ignore this fact at this point. Since 
sin 2 #14 is small, we do not expect big changes in the PDFs of |m ai g|, a, /3 = e, /i, r, except for the case \m ee \. We 
have explicitly checked that this is the case. 



In Fig. 10 we show \m ee \ and the new entries |m as |, a = e, /i, r for the normal hierarchy and Am 2 x = 1.71 eV 2 . 
For the case Am 2 x = 0.95 eV 2 , the plots would be similar in shape, however with different scale. For \m ee \ the largest 
value goes down from 0.06 eV to 0.045 eV. For |m es | the range changes from ~(0.12 - 0.26) eV to ~(0.09 - 0.20) eV. 
For Ira^l, \m TS \ there is basically no difference and for \m ss \ there is again a shift in the range from ~(1.25 - 1.30) 
eV to -(0.935 - 0.97) eV. 

In Fig. 11 we show the correlations among the PDF's of \m ee \ and |m es |, |m MS |, |m rs |, |m ss |, for the normal hierarchy 
and Am^~= 1.71 eV 2 . Since the difference between the case Am^ = 1.71 eV 2 and Am 2 x = 0.95 eV 2 is basically the 
scale, as commented above, we do not show the correlations in this case either. 

We observe that \m ee \ is very correlated with |m es |. This is easy to understand from the formulae in Appendix |A| 
as \m ee \ goes as s\ A while |m es | goes like 514 producing a squared root behavior. The thickness is driven by the CP 
phases. The linear behavior between \m ee \ and \m ss \ can be explained by noting that \m ss \ « 777,4 — \m ee \. Again due 
to the approximate \i-t symmetry we get similar ranges and behaviors for \m ee \ x \m^ s \ and \m ee \ x |m rs |. 

Also in this case it is interesting to consider future perspectives, in particular considering the interplay between 
oscillation experiments and cosmology. From the cosmology point a view, new data by Planck should be soon released, 
with a measurement of the number of light species. From the oscillation point of view, a part of the currently allowed 
parameter space (that allows for a solution to the reactor and Gallium anomalies) will be probed in the next couple 
of years by the NUCIFER [35] and STEREO [36] experiments. A much broader portion should be explored with a 
timescale of more than ten years by the updated version of KamLAND, CeLAND [37], in principle allowing to confirm 
or rule out the presence of a sterile neutrino independently from the cosmological measurements. 



V. FINAL DISCUSSION AND CONCLUSION 



In this paper we have analyzed in a probabilistic way the possible structures and correlations among the neutrino 
mass matrix elements in the standard neutrino oscillation framework, in view of the latest global analysis of the neu- 
trino oscillation data. This is done by constructing PDFs for each matrix element by assuming gaussian distributions 
for the known oscillation parameters and flat distributions for the unknown ones. 
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I Wee I (meV) \m efl \ (meV) 



FIG. 6: PDFs for the distribution of: \m ee \ x \m eT \ for the the normal hierarchy with sin 2 613 uncertainty reduced to a half of 
its present value (top left panel); \m^^\ x \m TT \ for the the normal hierarchy with Ara|i uncertainty reduced to a half of its 
present value (top right panel); \m ee \ x Ira^l for the the quasi-degenerate case with sin 2 612 uncertainty reduced to a half of 
its present value (bottom left panel) and \m e ^\ x |m MAt | for the the inverted hierarchy with sin 2 #23 uncertainty reduced to a 
half of its present value (bottom right panel). In all cases #23 was taken to be in the first octant. The colored areas are for the 
present uncertainties of the oscillation parameters, whereas the back lines are for the assumed future reduced uncertainty of 
one of the parameters. 

We analyzed the possible textures of the mass matrix allowed by data in the hierarchical and almost degenerate cases 
and discussed the future perspectives for better determining these matrix elements by future neutrino oscillation and 
non-oscillation data. The conclusion is that a better determination of the currently measured oscillation parameters 
will, in general, have a small effect on the matrix elements. The biggest effect will come from a better determination 
of sin 2 #23 solving the octant degeneracy, as one could have guessed. A determination of S would be significant, 
particularly for the normal hierarchy. In the inverted ordering, it would play a bigger role in the determination 
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FIG. 7: PDFs for the distribution of: \m ee \ x |m e/x | (top panels) and Ira^l x |m eT | (bottom panels) for the normal ordering 
From left to right 5 = 0°, 180° and 270°, assumed to be determined within 20°. 



of \m^ T \ and \m TT \. For the quasi-degenerate case, the impact of the determination of 5 is rather small due to the 
presence of more relevant CP phases. Future inputs from beta and neutrinoless double beta decay experiments, as well 
as cosmology, seem to be the most promising in providing new clues to understand the flavor structure. A specially 
encouraging scenario would be to have mp in the reach of KATRIN experiment, the neutrino mass hierarchy settled 
by near future oscillation experiments, and also a possible cosmological measurement of the sum of neutrino masses. 

Models of neutrino mixing based on discrete flavor symmetries, as discussed at the end of in Sec. |III[ anticipate 
relations among mixing angles and the 5 phase that can be tested in the future. If these relations turn out to be true, 
they will impose correlations among the mass matrix elements beyond the ones considered here. 

We extend our analysis to include the possibility of a sterile neutrino with mass and mixings allowed by the reactor 
and gallium anomalies (disregarding current cosmological bounds). We discuss what are the modifications to the mass 
matrix pattern in this scenario, finding only relevant modifications for m ee as expected. Despite the presence of more 
parameters than in the standard scenario, the larger sterile neutrino mass would make easier to measure both \m ee \ 
and mp. This may also have a big impact in cosmological models, as well as in the future strategies for oscillation 
experiments. 
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Appendix A: Matrix Elements Squared 



Here we give some approximate expressions for the matrix elements \m a p\ 2 . We use the following notation: x 2 = 
sin 2 #i2, y 2 = sin 2 #i 3 , z 2 = sin 2 623 and w 2 = sin 2 #14. Here we set yl — y 2 — » 1 and \/i — i/; 2 — » 1. The standard 
case is recovered by taking w —> 0. 



I |2 2 4, 2 4, 2/1 2\ 2 , 2 4 

|m ee | ~ m A w + m 2 x -\- m 1 [1 — x ) + m 3 y 

+ 2 mirri4 (l — x 2 ) iu 2 cos [2 (Ai — A4)] 

+ 2 m\rri2 x 2 (l — x 2 ) cos (2Ai) 

+ 2 mim 3 ?/ 2 (1 - x 2 ) cos [2 (6 - Ai + A 3 )] 

+ 2 m2m4 x 2 k; 2 cos (2A4) 

+ 2ra 2 ra 3 £ 2 2/ 2 cos[2(( s ) + A 3 )] 

+ 2 m 3 m4y 2 w 2 cos [2 (5 + A 3 - A 4 )] (AI) 
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FIG. 9: We show the current allowed regions for \m ee \ at 99% CL as a function of the effective electron neutrino mass, 777/3, 
on the left panel and as a function of the sum of the neutrino masses, ^rrii, on the right panel. The region allowed by the 
normal (inverted) mass ordering is in blue (red), the recent limit on \m ee \ given by KamLAND-Zen [29] in gray and the region 
excluded by cosmology [31] in magenta. We also show the reach expected for the beta decay experiment Katrin [32] , as well as 
the ultimate reach aimed by the neutrino less double beta decay experiments GERDA and CUORE according to Ref. [30] . 



I^e/il 2 ~ (l-x 2 )x 2 (l- z 2 )(m\ + m 2 2 ) + m\ (l - x 2 f y 2 + m\x^y 2 + m\y 2 
+ [2 m\m2 x 2 y 2 z 2 (l — x 2 ) — 2 m\m2 x 2 (l — z 2 ) (l — x 2 )] cos (2Ai) 
- 2m 1 m 3 y 2 z 2 (l-x 2 )cos[2(5- Ai + A 3 )] 

cos 5 



2 m\ xyz (l — x 2 ) 3 ^ 2 \J\ — z 2 — 2 m\ x 3 yz\/ 1 — x 2 y/l — z 2 

— 2 rri\rri2 xyz (l — x 2 ) 3 ^ 2 y/l — z 2 cos (6 — 2Ai) 

— 2 m2m3 x 2 y 2 z 2 COS [2 (5 + A3)] 

+ 2 m 2 7n 3 xyz\f\ — x 2 \/l — z 2 cos ((5 + 2A3) 

— 2 rn 1 rn 3 xyz\J\ — x 2 \/l — z 2 cos ((5 — 2Ai + 2A3) 
+ 2 mim2 x 3 yz\] 1 — x 2 \/l — z 2 cos (J + 2Ai) 



(A2) 



I 1 2 , / "i 2\ 2 2/ 2 1 2\ 1 2 A 2\ z 2 . 2 4 2, 2 2 /1 2\ 

|777 er I ~ (1 — X ) X Z (m 1 + 777 2 J + 7/7 x [1 — X J 7/ + 777 2 X 7/ + 7/7 3 7/ (1 - Z ) 

+ [2 777i7772 £ 2 7/ 2 (l — £ 2 ) (l — X 2 ) — 2 777i7772 X 2 ^ 2 (l — X 2 )] COS (2Ai) 

- 2 77717773 7/ 2 (l - Z 2 )(l - X 2 ) COS [2 (<5 - Ai + A3)] 

— 2 7/7 2 xt/z (l — x 2 ) 3 ^ 2 \/l — z 2 — 2 m\ x 3 yz\/ 1 — x 2 \/l — z 2 



cos J 



+ 2 77717772 X7/Z (l — X 2 ^^ \/ 1 — Z 2 COS ((5 — 2Ai) 

— 2 77727773 X 2 7/ 2 (1-Z 2 )C0S[2(5 + A 3 )] 

— 2 77727773 XyZ^J 1 — X 2 \/ 1 — Z 2 COS ((5 + 2A3) 

+ 2 rn 17773 xyzyj 1 — x 2 \/l — z 2 cos ((5 — 2Ai + 2A3) 

— 2 777i7772 x 3 yz\/ 1 — x 2 \/ 1 — Z 2 COS ((5 + 2Ai) 



(A3) 
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FIG. 10: PDFs for the distributions of the absolute value of |m ee |, \m e 
Amlx = 1.71 eV 2 . 



\m TS \ and |m ss |, for the normal hierarchy and 
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FIG. 11: PDFs for the distributions of the absolute values \m ee \ x |m es |, \Tn e e \ X Ira^sl, |m ee | x |m TS | and \m ee \ x |m ss |, for the 
normal hierarchy and Am^ = 1.71 eV 2 . 
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FIG. 12: We show the current allowed regions for \m ee \ at 99% CL as a function of the effective electron neutrino mass, m^, 
on the left panels and as a function of the sum of the neutrino masses, J^ m *? on ^ ne right panels for the 3+1 scenario. Here 
Aralx = 1.71 eV 2 . The region allowed by the normal (inverted) mass ordering is in blue (red), and the recent limit on \m ee \ 
given by KamLAND-Zen [29] in gray. We also show the reach expected for the beta decay experiment Katrin 32 , as well as 
the ultimate reach aimed by the neutrino less double beta decay experiments GERDA and CUORE according to Ref. [30] . 



\m^\ 2 w 4 (1 - z 2 ) x 2 y 2 z 2 (l - x 2 )(m\ + m%) + (1 - z 2 ) 2 [m\x A + m\(\ - x 2 ) 2 ] 

\2 



z 4 



2 (-\ 2\ z 4 i 2 4 4, 

if (1 - x z ) + m z 2 x*y* + 

+ 2 mim 2 x 2 (1 - x 2 ) [(1 - z 2 ) 2 - 4 ?/ V (l - z 2 ) + cos (2A X ) 

+ 2mim 2 (l-^ 2 )x 4 y 2 z 2 cos[2((5 + Ai)] 

+ 2 mim 2 (1 - z 2 )(l - x) 2 (l + x) Vz 2 cos [2 (5 - Ai)] 

+ 2 m 2 m 3 (1 - z 2 )(l - x 2 )z 2 cos (2A 3 ) 

+ 2 (m 2 + ml) (1 - z 2 )(l - x 2 )x 2 y 2 z 2 cos (25) 

+ 4 rrii'mz \f \ — z 2 \/l — x 2 xyz 3 cos (5 — 2Ai + 2A 3 ) 

— 4 m 2 m 3 \/l — \/l — x 2 xyz 3 cos (J + 2A 3 ) 

+ 4 mim 2 x 3 yz\/ 1 — z 2 \/l — x 2 [y 2 z 2 — (1 — z 2 )] cos (5 + 2Ai) 

- 4x^^/1 -^ 2 V 7 ! - ^ 2 [m? (x 2 ((y 2 + l) z 2 - l) - ?/V) 
+ ra 2 , (x 2 ((y 2 + 1) z 2 - 1) - ^ 2 + 1)] cos (S) 

+ 4 mim 2 E2/z<s/l - £ 2 \/l - ^ 2 (1 - x 2 ) [l - (1 + y 2 )z 2 ] cos (5 - 2Ai) 
+ 2 m 2 m 3 x 2 y 2 z A cos [2 ((5 + A3)] 

(1 - x 2 )?/V cos [2 (J - Ai + A 3 )] 
+ 2 mim 3 (1 - z 2 )x V cos [2 (Ai - A 3 )] (A4) 
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Im^l 2 « (l-z 2 )z 2 [m2 + (l- a; 2 ) 2 (m2y 4 + m2)+x 4 K + m2y 4 )] 

+ (1 - x 2 ) x 2 y 2 (m 2 + m 2 ) [(1 - z 2 ) 2 + z 4 } 

+ 2m im2 x 2 [z 2 {\ -z 2 ){{l + y 2 ) 2 - x 2 ) - y 2 {\ - x 2 )(z 4 + (1 - z 2 ) 2 )} cos (2A X ) 

+ 2 xyz(l — 2z 2 )yfl — z 2 \/l — x 2 [(1 — x 2 )(m 2 + m\y 2 ) — x 2 (m 2 + m\y 2 )\ cos 5 

- 2 m 1 m 2 xyz(l - 2z 2 )\/\ - z 2 (l - x 2 ) 3/2 (l + y 2 ) cos (S - 2A X ) 
+ 2mim 3 xyz(l - 2z 2 )\J\- z 2 \Jl - x 2 cos (6 - 2Ai + 2A 3 ) 

- 2 m 2 m 3 - 2z 2 )\A - z 2 y/l - x 2 cos (6 + 2A 3 ) 

- 2 {m\ + ml) (1 - z 2 )(l - x 2 )x 2 y 2 z 2 cos (25) 

- 2m 2 m 3 (l-z 2 )(l — x 2 )z 2 cos (2A3) 

+ 2 mim 2 x 3 ^(l - 2z 2 )y/l- z 2 y/l-x 2 (l + ?/ 2 ) cos (5 + 2Ai) 

- 2raira 2 2/V(l -z 2 )(l -x 2 ) 2 cos [2 (5- Ai)] 

- 2mim 2 xy^ 2 (l -z 2 ) cos [2 (5 + Ai)] 
+ 2m 2 m 3 a; 2 ^ 2 (l - z 2 ) cos [2 + A 3 )] 
+ mim 3 ?/ 2 (l - x 2 ) cos [2 (5 - Ai + A 3 )] 

- 2raira 3 (1 - z 2 )xVcos[2(Ai - A 3 )] (A5) 



4 (1 - z 2 ) x 2 y 2 z 2 (l - x 2 )(m{ + ra|) + z A [m\x A + m|(l - x 2 ) 2 ] 



+ 



mi 



2 /i 2\ 2 4 i 2 4 4, 2 /1 2\ 



^2\2 



+ 2mim 2 x 2 (1 -x 2 ) [y 4 (l - z 2 ) 2 -Ay 2 z 2 (l - z 2 ) + z 4 ] cos(2Ai) 

+ 2mim 2 (l-^ 2 )x 4 ?/ 2 z 2 cos[2(^ + Ai)] 

+ 2 mim 2 (1 - z 2 )(l - xf(l + x)V z 2 cos [2 (<J - Ai)] 

+ 2ra 2 ra 3 (l -z 2 )(l — x 2 )z 2 cos (2A 3 ) 

+ 2 (ra 2 + ra 2 ) (1 - z 2 )(l - x 2 )x 2 y 2 z 2 cos (26) 

— 4mim 3 \/l — z 2 \[\ — x 2 xyz(l — z 2 ) cos (5 — 2Ai + 2A 3 ) 
+ 4 m 2 m 3 \/l — ^ 2 \/l — x 2 xyz(l — z 2 ) cos (5 + 2A 3 ) 

+ 4mim 2 x 3 yz\/ 1 — z 2 \/l — x 2 [y 2 z 2 — y 2 + z 2 ] cos (5 + 2Ai) 

- Axyzy/l - x 2 \/l- z 2 [ra 2 ((l - x 2 ) ?/ 2 (l - z 2 ) + x 2 z 2 ) 



m 9 



(*V (l-z 2 )+z 2 (l-x 2 ))]cos(£) 



+ 4mim 2 -z 2 \J\- x 2 (1 - x 2 ) [(1 - z 2 )y 2 - z 2 ] cos (S - 2Ai) 

+ 2 ra 2 ra 3 x 2 t/ 2 (1 - z) 2 (l + z) 2 cos [2 (5 + A 3 )] 

+ 2 mxms (1 - x 2 )y 2 (l - z) 2 (l + z) 2 cos [2 (5 - A x + A 3 )] 

+ 2raira 3 (l -z 2 )x 2 z 2 cos[2(Ai - A 3 )] (A6) 



rajx 4 + ra 2 (l — x 2 ) + ra 2 ?/ 4 + ra 2 

— 2 ra 3 ra4 w 2 ?/ 2 cos [2 (5 + A 3 — A4)] 
+ 2 m\Vfi2 w 2 x 2 (l — x 2 ) cos (2Ai) 

+ 2 mim 3 w V (l - x 2 ) cos [2 (<S - Ai + A 3 )] 

— 2 mim 4 (l — x 2 ) cos [2 (Ai — A4)] 
+ 2 ra 2 ra 3 w 2 x 2 y 2 cos [2 (5 + A 3 )] 

— 2 m 2 m4 it; 2 x 2 cos (2A4) 



(A7) 
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\m s ^\ 2 « w 2 [ml (x 2 - 1) (x 2 ((y 2 + l) - l) - ?/V) + m 2 ^ 2 (x 2 ((y 2 + l) z 2 - l) - ^ 2 + l) + ra 3 2/V] 

+ 2 w 2 xyz\/ 1 — x 2 \/l — z 2 [(1 — x 2 )m\ — x 2 m^\ cos (5 

- 2mim 2 w;V(l-^) 3/2 V / l-^ cos (5-2Ai) 
+ 2 rri\<ni2 w 2 x 3 yz\/ 1 — x 2 \/l — z 2 cos (5 + 2Ai) 

— 2 m\m 3 w 2 xyz\/ 1 — x 2 \/l — z 2 cos ((5 — 2Ai + 2A3) 

+ 2 1712ms w 2 xyz\/ 1 — x 2 \/l — z 2 cos ((5 + 2A3) 

— 2 m^ms w 2 x 2 y 2 z 2 cos [2 (J + A3)] 

+ 2raira 2 wV(l -x 2 ) [(1 + y 2 )z 2 - l] cos (2Ai) 

- 2 mim 3 (1 - x 2 ) w 2 y 2 z 2 cos [2 (5 - Ai + A 3 )] (A8) 



|m sr | 2 « ^ 2 [m 2 (1 - x 2 ) ((1 - x 2 ) ^ 2 (1 - z 2 ) + xV) 

+ m 2 (xV + x 4 (^ 2 (1 - z 2 ) - z 2 )) + m 2 ^ 2 (l - z 2 )] 

— 2 w 2 xyz\/ 1 — x 2 \/l — z 2 [(1 — x 2 )m\ — x 2 m^\ cos 5 
+ 2 m\m2 w 2 xyz(l — x 2 ) 3 ^ 2 \/l — z 2 cos ((5 — 2Ai) 

— 2 m\m2 w 2 x 3 yz\/ 1 — x 2 \/l — z 2 cos ((5 + 2Ai) 

+ 2 rriirris w 2 xyz\/ 1 — x 2 \/l — z 2 cos (5 — 2Ai + 2A3) 

— 2 777,2 ra 3 w 2 xyz\/ 1 — x 2 \/l — z 2 cos ((5 + 2A3) 

- 2 m 2 m 3 w 2 x 2 y 2 (l - z 2 ) cos [2 (5 + A 3 )] 

- 2m 1 m 2 w 2 x 2 (l - x 2 ) [z 2 - (1 - z 2 )y 2 ] cos (2X 1 ) 

- 2m 1 m 3 (l-x 2 )w 2 y 2 (l-z 2 )cos[2(5-\ 1 +\ 3 )} (A9) 



|m ss | 2 ~ w A [(1 — x 2 ) 2 m\ + x 4 m?, + y 4 m 2 ] + m 2 

+ 2 m\m2 w A x 2 (l — x 2 ) cos(2Ai) 

+ 2 77727774 w 4 x 2 cos(2A4) 

+ 2m 1 m 3 w A y 2 (l - x 2 ) cos [2 (S - X 1 + A 3 )] 

+ 2 777,1777,4 w 2 (l — x 2 ) cos [2 (Ai — A4)] 

+ 2 777,2777,3 w 4 x 2 y 2 cos [2 (J + A3)] 

+ 2 m 3 m 4 w 2 y 2 cos [2 (J + A 3 - A 4 )] (A10) 
Appendix B: Complete Set of Correlation Plots for the Matrix Elements 
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FIG. 13: PDFs for the distribution of the absolute value of several pairs of matrix elements. We use blue, green and red for 
the allowed region at 68.27%, 95.45% and 99.73% CL, respectively. Here mi —> and 623 is assumed to be in the first octant. 
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FIG. 14: Same as Fig[l3]but for #23 in the second octant. 
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FIG. 15: PDFs for the distribution of the absolute value of several pairs of matrix elements. We use blue, green and red for 
the allowed region at 68.27%, 95.45% and 99.73% CL, respectively. Here — »• and #23 is assumed to be in the first octant. 



23 



30 

> 

CD 

£20 
10 



' i 1 1 1 1 i 1 




20 30 40 50 
|m ee | (meV) 



30 
25 
| 20 
— 15 

5 



i ■ 1 ■ 1 1 1 ■ 1 ■ i ■ 1 ■ 1 1 1 




10 20 30 40 50 
|m ee | (meV) 




10 20 30 40 50 
|m ee | (meV) 



40 
>30 

CD 

^20 

_|io 





' I 1 1 1 1 I 1 




20 30 40 50 
|m ee | (meV) 



10 20 30 40 50 
|m ee | (meV) 



30 



& 20 



A 10 



0- 



■ I 1 1 1 1 I 1 




40 



^30 
> 

CD 

^20 
-^10 




10 20 30 40 
|m e/ J (meV) 



30 
25 

I 20 

— 15 

i,o 

5 



'i 1 1 ■ ■ i 1 ■ 1 ■ i 1 1 ■ ■ i ■ 




10 20 30 
\m efi \ (meV) 



10 20 30 
\m efl \ (meV) 



40 
>30 

CD 

£20 

s 10 




10 20 30 
|m e/ J (meV) 




30 

CD 

B 20 



10 




10 20 
\m m \ (meV) 



10 20 30 40 
|m eT | (meV) 





10 20 30 40 
|m eT | (meV) 




10 20 
\m m \ (meV) 



1 20 

— 15 




10 20 
\m m \ (meV) 




5 10 15 20 25 30 
\m^\ (meV) 



FIG. 16: Same as Fig[l5]but for 623 in the second octant. 




FIG. 17: PDFs for the distribution of the absolute value of several pairs of matrix elements. We use blue, green and red for 
the allowed region at 68.27%, 95.45% and 99.73% CL, respectively. Here mi =0.1 eV, we impose the normal mass ordering 
and 623 is assumed to be in the first octant. 
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FIG. 18: Same as Fig[l7|but for 623 in the second octant. 



